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In the standard model, point like parton-parton scatterings in high energy proton-proton colli- 
sions can give rise to final states with energetic jets. At large momentum transfers, events with 
at least two energetic jets (dijets) may be used to confront the predictions of perturbative Quan- 
tum Chromodynamics (pQCD) and to search for signatures of new physics. In parton-parton 
scattering, the angular distribution of the outgoing partons, da /d cos 9*, is directly sensitive 
to the spin of the exchanged particle, where a is the parton-level cross section and 9* is the 
polar scattering angle in the parton-parton center-of-mass (CM) frame. While QCD predicts 
a noticeable deviation of the dijet angular distribution from Rutherford scattering, at small 
CM scattering angles the angular distribution is proportional to the Rutherford cross section, 
da/dcos 9* ~ 1/(1 — cos 9*) 1 , characteristic of spin-1 particle exchange. The dijet angular dis- 
tributions do not strongly depend on the details of the parton distribution functions (PDFs), 
since the angular distributions for the underlying processes, qg —> qg, qq' — > qq' , and gg — > gg, 
are similar. 

For the scattering of massless partons, which are assumed to be collinear with the beam pro- 
tons, the longitudinal boost of the parton-parton CM frame with respect to the proton-proton 
CM frame, J/boost/ an d 9* are obtained from the rapidities \j\ and 1/2 of the jets from the two 
scattered partons by i/boost = 5(1/1 + Vi) and | cos0*| = tanhy*, where y* = \\y\—yi\ and 
where ±y* are the rapidities of the two jets in the parton-parton CM frame. The rapidity 
is related to the jet energy E and the projection of the jet momentum on the beam axis p z 
by y = jln[(£ + p z )/ (E — p z )]. The variable Xdijet = ex p(2y*) is used to measure the di- 
jet angular distribution, which for collinear massless-parton scattering takes the form Xdijet = 
(1 + I cos (1 — I cos 0*|). This choice of Xdi\et, rather than 9* , is motivated by the fact that 
^Cdijet/^dijet is flat for Rutherford scattering. It also allows signatures of new physics that 
might have a more isotropic angular distribution than QCD (e.g. quark compositeness) to be 
more easily examined as they would produce an excess at low values of ^dijet- The quantity 
studied in this analysis is (1 / Cdijet) (^dijet/ ^/tdijet)/ f° r several ranges of the dijet invariant mass 
Mjj. Previous searches for quark compositeness using the dijet angular distribution or related 
observables in pp and pp collisions have been reported at the SppS by the UA1 [lj collabo- 
ration, at the Tevatron by the DO HI [3J and CDF [4] collaborations, and at the Large Hadron 
Collider (LHC) by the ATLAS [5] collaboration. The CMS collaboration has also published a 
search on quark compositeness with a smaller data sample using the dijet centrality ratio 0. 
In this Letter, we present the first measurement of dijet angular distributions from CMS in pp 
collisions at y/s = 7 TeV. 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing an axial field of 3.8 T. Within the field volume are the silicon pixel and silicon 
strip tracker, the electromagnetic calorimeter (ECAL) and the hadron calorimeter (HCAL). The 
ECAL is made up of lead-tungstate crystals, while the HCAL is made of layers of plates of 
brass and plastic scintillator. These calorimeters provide coverage in pseudorapidity up to 
\t] I < 3, where rj = — lntan(0/2) and 9 is the polar angle relative to the counterclockwise 
proton beam direction. An iron/ quartz-fiber Cerenkov hadron calorimeter (HF) covers pseu- 
dorapidities 3 < \t]\ < 5. In addition, a preshower detector made of silicon sensor planes 
and lead absorbers is located in front of the ECAL at 1.653 < \tj\ < 2.6. The calorimeter 
cells are grouped in projective towers of granularity in pseudorapidity and azimuthal angle of 
0.087 x 0.087 at central pseudorapidities, with coarser granularity at forward pseudorapidities. 
Muons are measured in gas-ionization detectors embedded in the steel magnetic field return 
yoke. A detailed description of the CMS detector can be found elsewhere 0- 

Events were collected online with a two-tiered trigger system: Level-1 (LI) and the High Level 
Trigger (HLT). For this study, events were selected with five inclusive single-jet triggers, with 
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the following jet transverse momentum pj thresholds at LI (HLT): 20 GeV (30 GeV), 30 GeV 
(50 GeV), 40 GeV (70 GeV), 60 GeV (100 GeV), and 60 GeV (140 GeV). The jets at LI and HLT 
were reconstructed using energies measured by the ECAL, HCAL, and HF, and were not cor- 
rected for the jet energy response of the calorimeters. All except the highest-threshold jet trigger 
were prescaled as the LHC instantaneous luminosity increased during the course of data tak- 
ing. In each case, the trigger efficiency was measured as a function of dijet invariant mass Mo- 
using events selected by a lower-threshold trigger. For the analysis, Mjj- and ^dijet" re gi° ns were 
chosen such that the trigger efficiencies exceeded 99%. 

Jets were reconstructed offline from energies measured in the calorimeter towers using the anti- 
kj clustering-algorithm [8] with a distance parameter R = 0.5. Spurious jets from noise and 
non-collision backgrounds were eliminated by loose quality criteria on the jet properties 0. 
The jet four-momenta were corrected for the non-linear response of the calorimeters [10]. The 
performance of the CMS detector with respect to jet reconstruction is described in detail else- 
where ITU . 

Events were required to have a primary vertex reconstructed within 24 cm of the detector center 
along the beam line [12]. Events having at least two jets were selected and the two highest- 
Pt jets were used to measure the dijet angular distributions for different ranges in Mjj. We 
required ^dijet < 16 and |j/boostl < 1-H/ thus restricting the rapidities j/i and 1/2 of the two 
highest-pj jets to be less than 2.5. Nine analysis ranges were defined with the boundaries 
0.25 < Mjj < 0.35 TeV, 0.35 < Mjj < 0.5 TeV, 0.5 < Mjj < 0.65 TeV, 0.65 < Mjj < 0.85 TeV, 
0.85 < Mjj < 1.1 TeV, 1.1 < Mjj < 1.4 TeV, 1.4 < Mjj < 1.8 TeV, 1.8 < Mjj < 2.2 TeV, and 
Mjj > 2.2 TeV. The data correspond to integrated luminosities of 0.4, 3.5, 9.2, and 19.8 pb _1 
for the lowest four Mjj ranges and 36 pb _1 for the remaining ones. The uncertainty on the 
integrated luminosity has been estimated to be 11% fll3fl. 

The dijet angular distributions are corrected for migration effects in ^dijet and Mjj due to the fi- 
nite jet energy and position resolutions of the detector. The correction factors were determined 
using two independent Monte Carlo (MC) samples: pythia 6.422 HH with tune D6T Jl5]| and 
HERWIG++ 2.4.2 [16]. The four-momentum, rapidity, and azimuthal angle of each generated jet 
were smeared to reproduce the measured resolutions. The ratio of the two dijet angular distri- 
butions (the generated distribution and the smeared one) determined the unfolding correction 
factors for a given MC sample and for each Mjj range. The average of the correction factors for 
each Mjj range from the two MC samples formed the final unfolding correction applied to the 
data. The correction factors change the normalized dijet angular distributions for all Mjj ranges 
by less than 3%. For each Mjj range, the systematic uncertainty associated with each correction 
factor was set at 50% of its value. This approach covers the variations of the unfolding correc- 
tion factors determined from HERWIG++ and different PYTHIA tunes (D6T and Z2 [17]) that 
vary on their modelling of the jet kinematic distributions. The use of a parameterized model 
to simulate the finite jet pj and position resolutions of the detector, to determine the unfolding 
correction factors, resulted in a systematic uncertainty. This was estimated to be less than 1% 
for all Mjj ranges and was added in quadrature to the unfolding uncertainties. 

The normalized dijet angular distributions are relatively insensitive to many systematic effects, 
in particular they show little dependence on the overall jet energy scale. However, since Xdi]et 
depends on y* , they are sensitive to the rapidity dependence of the jet energy calibration. Typ- 
ical values for the jet energy scale uncertainties for the considered phase space in the variables 
of jet pr and rj covered in this analysis are between 3% and 4% fllOl . The uncertainty on the ^dijet 
distributions due to the jet energy calibration uncertainties was found to be less than 2.5%. The 
uncertainty on the dijet angular distributions from the jet pj resolution uncertainty, estimated 
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to be 10% IfTTl , was found to be less than 1%. The total systematic uncertainty on the Xdijet 
distributions, calculated as the quadratic sum of the contributions due to the uncertainties in 
the jet energy calibration, the jet pj resolution, and the unfolding correction, is less than 3% for 
all Mjj ranges. 

The corrected differential dijet angular distributions for different M« ranges, normalized to 
their respective integrals, are shown in Fig. [T] The data are compared to pQCD predictions at 
next-to-leading order (NLO) calculated with NLOJET++ |18| in the FASTNLO |19] framework. 
The calculations were performed with the CTEQ6.6 PDFs 1 20 1 . The factorization (}i A and renor- 
malization {y. r ) scales were set to (pr), the average dijet pj. Non-perturbative corrections due 
to hadronization and multiple parton interactions, determined using the average correction 
from PYTHIA (D6T tune) and HERWIG++, were applied to the prediction. The uncertainties on 
the pQCD predictions, indicated by the shaded band in Fig. [TJ are less than 6% (9%) at low 
(high) Mjj. These uncertainties include contributions due to scale variations and PDF uncer- 
tainties, as well as the uncertainties from the non-perturbative corrections. The uncertainty due 
to the choice of }if and }i r scales was evaluated by varying the default choice of scales in the 
following six combinations: (p/, p r ) = ((pr)/2, (pr)/2), ((pr)/2, (Pt)), {{Pt), (Pt)/2), (2(pr), 
2(Pt))/ (2(pr)/ (Pt)) an d ((pr)/ 2(pr))- These scale variations modify the predictions of the 
normalized Xdijet distributions by less than 5% (9%) at low (high) Mjj. The uncertainty due to 
the choice of PDFs was determined from the 22 CTEQ6.6 uncertainty eigenvectors using the 
procedure described in Ref. fl20"|, and was found to be less than 0.5% for all Mu ranges. Half the 
difference between the non-perturbative corrections from PYTHIA and HERWIG++ was taken as 
the systematic uncertainty, and was found to be less than 4% (0.1%) at low (high) Mu. Overall 
there is good agreement between the measured dijet angular distributions and the theoretical 
predictions for all Mu ranges. 

The measured dijet angular distributions can be used to set limits on quark compositeness 
represented by a four-fermion contact interaction term in addition to the QCD Lagrangian. 
The value of the mass scale A characterizes the strengths of the quark substructure binding 
interactions and the physical size of the composite states. A color- and isospin-singlet con- 
tact interaction (CI) of left-handed quarks gives rise to an effective Lagrangian term: C qq = 
?/o(27r/A 2 )(^ L 7?'^L)(^£7fi^L) fl2~l~ll22| , where rjo = +1 corresponds to destructive interference 
between the QCD and the new physics term, and tjo = —1 to constructive interference. We in- 
vestigate a model in which all quarks are considered composite as implemented in the PYTHIA 
event generator. 

The contributions of the CI term in PYTHIA are calculated to leading order (LO), whereas the 
QCD predictions for the dijet angular distributions are known up to NLO. In order to account 
for this difference in the QCD plus CI prediction, the cross-section difference ~~ <7 wi D 

was added to the LO QCD+CI prediction in each Mjj and Xdijet bin. With this procedure, we 
obtain a QCD+CI prediction where the QCD terms are corrected to NLO while the CI terms 
are calculated at LO. Non-perturbative corrections due to hadronization and multiple parton 
interactions were also applied to the prediction. The prediction for QCD+CI at the scale of 
A + = 5 TeV (t]o = +1) and A~ = 5 TeV (rjo = — 1) are shown in Fig. [TJ for the four highest Mjj 
ranges. 

We perform a statistical test discriminating between the QCD-only hypothesis and the QCD+CI 
hypothesis as a function of the scale A based on the log-likelihood-ratio Q = — 21n( L ^°^ CI ). 
The likelihood functions Lqcd+ci an d Lqcd are modelled as a product of Poisson likelihood 
functions for each bin in ^dijet and Mjj in the four highest Mjj ranges. The prediction for each 
Mjj range is normalized to the number of data events in that range. The p-values, Pqcd+ci(Q > 
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Q bs) and Pqcd(Q < Qobs)/ are obtained from ensembles of pseudo-experiments. A modified 
frequentist approach p3ti25l based on the quantity 

p L PQCD+Cl(Q > Qobs) 
S ~ 1 ~ PqCv(Q < Qobs) 

is used to set limits on A. This approach is more conservative than a pure frequentist approach 
(Neyman construction) and prevents an exclusion claim when the data may have little sensitiv- 
ity to new physics |26|. Systematic uncertainties were introduced via Bayesian integration ||27| 
by varying them as nuisance parameters in the ensembles of pseudo-experiments according to 
a Gaussian distribution convoluted with the shape variation induced to the ^dijet distributions. 
We consider the QCD+CI model to be excluded at the 95% confidence level if CL S < 0.05. Fig- 
ure^ shows the observed and expected CL S as a function of the CI scale A + . From this we 
determine the lower limit on A + to be 5.6 TeV. The observed limit agrees within 1.4 standard 
deviations with the expected limit of 5.0 TeV, which was evaluated at the median of the test 
statistics distribution of the QCD model. The observed limit is slightly higher than the ex- 
pected one because, for the range M» > 2.2 TeV, the measured dijet angular distribution at low 
Xdijet i s lower than, although statistically compatible with, the QCD prediction. The limit for 
the CI scale was also extracted using an alternate procedure in which the data were not cor- 
rected for detector effects and instead the MC predictions were resolution-smeared. The limit 
obtained was found to agree with the quoted one within 0.4%. The corresponding observed 
and expected limits on A~ are 6.7 TeV and 5.8 TeV, respectively. 

Shortly before the completion of this Letter, an exact NLO calculation of QCD effects to quark 
compositeness became available [28 j. This calculation indicates that the limit on A + obtained 
in the present analysis, which only takes into account the LO prediction for the contribution of 
the contact interaction, might be overestimated by up to 10% compared to the value obtained 
if the NLO calculation were used. 

In summary, CMS has measured the dijet angular distributions over a wide range of dijet in- 
variant masses. The Xdijet distributions are found to be in good agreement with NLO pQCD 
predictions, and are used to exclude a range of a color- and isospin-singlet contact interaction 
scale A for a left-handed quark compositeness model. With a modified frequentist approach, 
a lower limit on the contact interaction scale of A + = 5.6 TeV (A - = 6.7 TeV) for destructive 
(constructive) interference at the 95% confidence level is obtained, which may be compared 
with a limit of 5.0 TeV (5.8 TeV) expected for the number of events recorded. These are the 
most stringent limits on the contact interaction scale of left-handed quarks to date. 
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Figure 1: Normalized dijet angular distributions in several Mjj ranges, shifted vertically by the 
additive amounts given in parentheses in the figure for clarity. The data points include statis- 
tical and systematic uncertainties. The results are compared with the predictions of pQCD at 
NLO (shaded band) and with the predictions including a contact interaction term of compos- 
iteness scale A + = 5 TeV (dashed histogram) and A~ = 5 TeV (dotted histogram). The shaded 
band shows the effect on the NLO pQCD predictions due to y. r and \ij scale variations and 
PDF uncertainties, as well as the uncertainties from the non-perturbative corrections added in 
quadrature. 
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A + [TeV] 

Figure 2: Observed CL S (solid line) and expected CL S (dashed line) with one (two) standard 
deviation(s) indicated by the dark (light) band as a function of the contact interaction scale A + . 
The 95% confidence level limits on A + are extracted from the intersections of the observed and 
expected CL S lines with the horizontal line at CL S =0.05. 
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